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SUMMARY 


The basic whirl stability of a rotating windmill on a flexible tower is 
reviewed. Effects of unbalance, gravity force, gyroscopic moments, and aero- 
lynamics are discussed. Some experimental results on a small model windmill 
ire given. 

INTRODUCTION 


There has been a renewed interest in the use of large windmills for 
generating power. Such large, rotating structures mounted on tall flexible 
rowers may give rise to significant vibration and fatigue problems. A good 
leal of the experience and knowledge gained during the last few years in con- 
nection with helicopter rotors and tilt-wing proprotors can be applied to such 
Large windmill systems. However, there are unique features of windmills and 
their operating environment that will have to be explored individually. 

A basic description of general rotating machinery problems can be found in 
Den Hartog’s book, (ref. 1). Loewy (ref. 2) presents a good review of rotary 
tfing dynamic and aeroelastic problems. More recently, a NASA special publica- 
tion (ref. 3) gives a good sampling of current problems dealing with rotor 
dynamics. References 4, 5, 6 are typical of recent investigations of problems 
tf large windmill systems. The present article will first review some dynamic 
problems of a rotating windmill on a flexible tower, then present some pre- 
liminary experimental results on a small windmill model. 


REVIEW OF THEORY 


Figure 1 shows the model used *for representing a windmill rotor mounted on 
a flexible tower. There is an absolute axis system x, y, z fixed in space, 
and also an axis system x s , ys, z s along the windmill shaft and having x s lie 
Ln the vertical plane (plane of xz) . The ith blade rotates about the axis z s 
tfith a constant speed fl, and can lag an angle <t>± in x s y s , plane and flap an angle 
3-j_ perpendicular to x s y s plane. Any point, on the blade can be expressed 
relative to the shaft axes x g , y g , z g as 
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x = e cos il) . + E, cos OK + 4.) cos 3. 

S r l r l T x 1 

y = e sin ll) . + E, sin OK + <b .) cos 3. 

J’g T x ’ r x r x x 

z g = ? sin 3 

In the above, represents the angular position of the i^ blade and e is, the 
hinge off-set. The origin of the shaft axis is assumed to translate fore-and- 
aft a distance qp and laterally a distance qp. Associated with these deflec- 
tions are an angular rotation 0pqp about the y s axis, another possible rota- 
tion Spqp about the x s axis, and a vertical deflection h v qp in the x direction. 
The coefficients 9p, 0p, by can be obtained from the vibration modes of the 
tower (often, h v ~ -h p) . The shaft axes can. be located relative to the fixed 
axes by performing a rigid body rotation about the y s axis and about the x g 
axis respectively. This gives the relation 
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Using the small angle approximation, sin 0^qp ^ 0^qp, cos 0pqp * 1 - 0 q /2 etc, 
in equation (2) and combining with equation (1) and the appropriate deflections 
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where x g , y s , z g are given by equation (1) . The velocity components x, y, z ■ 
are obtained from equations (3) '.by differentiation with respect to time t. 

Then, by forming the kinetic energy of the blades and tower, and placing into 
Lagrange’s equations, one can obtain the equations of motion of the windmill 
system. To simplify the lengthy algebra involved, it was assumed the hinge 
offset e = 0, and only those terms leading to linear terms in the final equa- 
tions of motion were retained. The following standard mass integrals were de- 
fined for the i 1 -* 1 blade, 

M i = J dm ’ S i = dm » J ± = J^ 2 dm (4) 

In the development, a two-bladed rotor was assumed with slightly unequal 
masses, such that Mp = Mg + M u /2 and M 2 = Mg - M u /2 where Mg was the average 
mass and My the unbalance in mass of the blades. Similar definitions were 
made for the average and unbalance in moment Sg and S u , and in moment of 
inertia Ig and I u . The vertical gravity loads were put in by writing the 
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incremental work as, 


6 W = J [f x <5x + f 5y + f z 6y]d£ = ^n^^n (5) 

where f x = -mg, f = f z = 0 , 6q n represents <Sq p , 6qx , <$3i,<S<J>i respectively, 
and 6x, 6y, Sz are fbund by differentiating equation (3). A similar procedure 
could be used for obtaining the aerodynamic forces acting on the blade. How- 
ever there, it is convenient to relate the air forces perpendicular and paral- 
lel to the blade axis §, 

The final, linear equations of motion in terms of the six coordinates q F , 
<lL» $1> $2> 4»i* $2 are » 
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* 

n the above equations, the k n q n and c n q n terms represent structural stiffness 
nd damping, the g terms represent the effect of gravity loads, and the Q 
erms represent the aerodynamic forces. The My F and are the generalized 
ower masses corresponding to q p and q L respectively. 
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Some of the gravity loads act as stiffness terms in the equations. The blade 
coordinates i|^ = fit and ^2 “ $1 + Tr * For t ^ ie two-bladed case, it is sometimes 
convenient to introduce the symmetric and antisymmetric blade variables, 

e s - + 3 2 )/2> 3 a “ (3 i " 3 2 )/2> V = etc * do.) 

to lessen the coupling between the degrees of freedom. Indeed, for a com- 
pletely balanced rotor without gravity effects, the 4> s would be uncoupled from 
the other equations. In general though, all six coordinates are involved. 

Equations (6) to (9) are a linear set of equations with periodic co- 
efficients, subjected to gravity, rotor unbalance S u , and aerodynamic wind 
forcing functions. The gravity loads act directly on the blades while the un- 
balance loads shake the tower which in turn couples into the blades. In addi- 
tion to forced response, the homogeneous equations themselves may have strong 
instabilities present. These are generally investigated by the use of Floquet 
theory for these periodic coefficient equations. It should also be mentioned 
that for a three or more bladed rotor, the analysis is generally easier since 
one can eliminate the periodic coefficients by a suitable transformation of 
coordinates (at least for the balanced rotor, without gravity effects). See 
for example reference 7. 

Various investigators have examined different subcases of equations (6) to 
(9). Coleman and Feingold (ref. 8) first looked at the case qp = 0, 0^ = 0, 

0L = 0, with no gravity, unbalance, or aerodynamics present. Strong mechanical 
instabilities of a whirling nature were found to be possible at certain rota- 
tional speeds ■, involving coupling of lateral motion q^ with lag angle <j>^. 

This is the so-called "ground resonance" helicopter phenomenon. Reed (ref. 9) 
looked at the case 0^ = 0, <j>i - 0 with aerodynamics present. Again, strong 
instabilities were found involving q^ and the vertical hyqp coupling through 
the mechanical and aerodynamic gyroscopic forces (fiqp, fiqj, terms). This is the 
so-called "propellor whirl" flutter. Young and Lytwyn (ref. 10) looked at the 
case <j>^ = 0 with aerodynamics present. This is essentially "propeller whirl" 
with flapping. Johnson (ref. 11) has looked in detail at the whole coupled 
system, but without gravity and unbalance effects in connection with his studie 
of proprotors. Equations very similar to the ones here are presented there. 
Finally, it should be mentioned there is a whole series of detailed investiga- 
tions of rotors attached to fixed hubs (qp = 0, q^ = 0) which emphasize the 
aerodynamic interaction between blade flapping, lagging, pitching and nonlinear 
dynamic effects brought on by large initial coning angles for the blades. See 
for example, references 4, 5, and 6. 


EXPERIMENT 


Some preliminary tests were run on a small .915 m (3.0 ft) diameter wind- 
mill placed in a wind tunnel. The general layout is shown in figure 2. The 
windmill had generally 2 blades , cantilevered in both the flap and lag direc- 
tions. The approximately uniform, untwisted blades had a .0762 m (3 in) chord, 
and could be set at any incidence angle. For a few runs, 4 blades were 
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attached to the windmill. 

The weight of a typical blade was .175 kg (.386 lbs). The cantilever 
natural frequencies of the non-rotating blades were measured as 33, 93, 172, 
and 310 Hz for the 1 st flap bending, 1 st lag bending, 2 nc ^ flap bending, and 
1 st torsion modes respectively. These were corrected for rotational effects in 
the standard manner, 0)^ = oo|^ + Lf2 2 , to give the rotating natural frequencies 
shown in figure 3. The tower stand had natural frequencies of 8.8, 16, 25 and 
75 Hz for the lateral yawing, vertical pitching, lateral translation and verti^ 
cal translation modes respectively. The windmill was instrumented to measure 
flap and lag bending moments at the blade root, and also lateral and vertical 
accelerations of the tower near the front bearing. 

The wind tunnel was run to about 18 m/sec (59.1 ft/sec), and after taking 
iata on windmill performance, the wind was turned off and the windmill would 
coast down to zero rotational speed. This gave a continuous frequency record 
Lhrough all the resonances of the system. Figures 4, 5, and 6 show the 
neasured bending moments and accelerations from such sweeps for a blade setting 
angle 0 =0°. Many superharmonic resonances can beTseen for the flap and lag 
sending moments. These occur near integer orders of the rotation frequency as 
zan be seen from figure 3. Particularly strong vibrations occured at 2 per 
revolution for both flap and lag. Indeed, lag moments near 10 times the static 
gravity moments are present at 50 Hz. The corresponding accelerations show a 
strong lateral resonance near 24 Hz. In these tests there was a small static 
inbalance due to unequal blade weights. Subsequent tests with another set of 
blades having a greater unbalance showed the same vibration patterns, but 
tfith peak amplitudes increased more than double. Also, tests run with four 
blades on the rotor showed similar strong resonances at 2 per revolution. The 
strong resonances in figures 4 to 6 seem then to have been caused by the 
rotating unbalance of the blade exciting tower stand frequencies which in turn 
excite blade frequencies superharmonically . Further details of these tests 
zan be found in reference 12. 


CONCLUSIONS 


A brief review of some of the dynamic problems associated with large 
rotating windmills has been given, together with some preliminary experimental 
results. The importance of flexible towers and their interaction with the 
rotating blade dynamics has been discussed. Although much work has already 
>een done in this area, many interesting dynamic problems remain to be re- 
solved, particularly those Involving rotors with built-in coning angles where 
xonlinear dynamics must be considered. 
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Figure 3.- Rotating natural frequencies of blades* 



Figure 4.- Flap bending moment vibrations. 





